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ABSTRACT 

We continue our investigation of overcoming Gibbs phenomenon, i.e., to obtain exponen- 
tial accuracy at all points (including at the discontinuities themselves), from the knowledge 
of a spectral partial sum of a discontinuous but piecewise analytic function. We show that if 
we are given the first N Gegenbauer expansion coefficients, based on the Gegenbauer polyno- 
mials C£{x) with the weight function (1 - x 2 )^* for any constant g > 0, of an L i function 
f(x), we can construct an exponentially convergent approximation to the point values of 
f(x) in any sub-interval in which the function is analytic. The proof covers the cases of 
Chebysliev or Legendre partial sums, which are most common in applications. 


1 Research supported by AFOSR grant 93-0090, ARO grant DAAL03-91-G-0123, DARPA grant N00014- 
91-J-4016, NSF grant DMS-921 1820, NASA grant NAC.1-1145, and NASA contract NAS1-19480 while the 
authors were in residence at ICASE, NASA Langley Research Center, Hampton, VA 23681. 
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1 Introduction 


In this paper we continue our investigation of overcoming the Gibbs phenomenon, i.e., recov- 
ering pointwise exponential accuracy at all poirits including at the discontinuities themselves, 
from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic func- 
tion, which we started in [4], [5], and [6]. 

Spectral approximations, such as the Fourier approximation based upon trigonometric 
polynomials for periodic problems, and the Chebyshev, Legendre or the general Gegenbauer 
approximation based upon polynomials for non-periodic problems, are exponentially accurate 
for analytic functions [3], [2]. However, for discontinuous but piecewise analytic functions, 
the spectral partial sum approximates the function poorly throughout the domain. Away 
from the discontinuity only first order accuracy is achieved. Near the discontinuity there are 
0(1) oscillations which do not decrease with N, the number of terms retained in the spectral 

sum. This is known as the Gibbs phenomenon. 

Our framework in [4], [5] and [6] to overcome Gibbs phenomenon to obtain exponential 
accuracy at all points for piecewise analytic functions relies heavily on using the Gegenbauer 
polynomials C x (x), which are orthogonal in [—1,1] with the weight function (1 — x 2 ) x for 
large A. We assume that the first -N < k < N Fourier coefficients, or the first 0 < k < N 
Legendre coefficients, of a discontinuous but piecewise analytic function, are given. The 
procedure consists of two steps: 

1. Using the given spectral partial sum of to recover the first m ~ N Gegenbauer 
expansion coefficients, based on a sub-interval [a, 6] C [—1,1] in which the function 
is presumably analytic, with exponential accuracy. This can be achieved for 
any L\ function, as long as we choose A in the weight function of Gegenbauer 
polynomials to be proportional to N. The error incurred at this stage is called 
the truncation error. 

2. For an analytic function in [a, b], proving the exponential convergence of its Gegen- 
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bauer expansion, when the parameter A in the weight function is proportional to 
the number of terms retained in the expansion. The error at this stage is labeled 
the regularization error . 


In [6] we demonstrated this procedure in the case of a discontinuous but piecewise analytic 
function, provided its Fourier or Legendre spectral partial sum is given. 

The proof of the Legendre case in [6] is based upon first expanding the Legendre poly- 

i_ 

iiomial Pk{x) = C£(x) into its Fourier series. It was essential in this proof that the Fourier 
expansion for the Legendre polynomial Pn(x), for large N, contains lower terms that decay 
exponentially with N (formula (2.13) in [6]). Unfortunately, it seems that this fact is true 
only for Legendre polynomials, probably because their weight function is special (= 1). It 
seems not true for other Gegenbauer polynomials, such as G'hebyshev polynomials. In an 
earlier version of [6], we quoted a formula (7.354, page 836 of [7]) to this effect for Ohebyshev 
polynomials. However, it is doubtful that Formula 7.354 of [7] is correct. 

In this paper, we will consider the case of general Gegenbauer spectral methods, with 
G’hebyshev and Legendre methods as special cases. We assume that f(x) is an Lj function 
on [-1,1] and analytic in a subinterval [a, 6] C [-1, 1]. We also assume that the Gegenbauer 
partial sum of f(x), based upon the Gegenbauer polynomials C£{x) with the weight function 
(1 — x 2 ) ,i 2 for any constant /x > 0, over the full interval [—1, 1], is known. The objective is 
to recover exponentially accurate point values over the subinterval [a, b] of analyticity. 

We will follow the same path as in [6]. Basically we will show that the first 0 < k < N 
Gegenbauer expansion coefficients, based on the Gegenbauer polynomials C£(x) for any 
constant /i > 0, contain enough information, such that a different, rapidly converging 
Gegenbauer expansion in the subinterval [a, 6], with the parameter A in the weight func- 
tion (1 — £ 2 ) A- * being proportional to N, can be constructed. As before, we will separate 
the analysis of the error into two parts: truncation error and regularization error. Trunca- 
tion error measures the difference between the exact Gegenbauer coefficients with A ~ N, 
and those obtained by using the spectral partial sum. This will be investigated in Section 3. 
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The regularization error measures the difference between the Gegenbauer expansion using 
the first few Gegenbauer coefficients with A ~ N, and the function itself in a sub-interval 
[a,b\, in which the function is assumed analytic. This error is estimated in [6] and we wdl 
simply quote the result in Section 4. The results are summarized in Theorem 4.3, and some 
remarks are also given in Section 4. Section 5 contains two numerical examples to illustrate 
our results. In Section 2 we collect some useful properties of Gegenbauer polynomials to be 
used later. 

Throughout this paper, we will use A to denote a generic constant or at most a polynomial 
in the growing parameters, as will be indicated in the text. It may not be the same at different 
locations. 


2 Preliminaries 


In this section we collect some useful results about the Gegenbauer polynomials, to be used 
in later sections. We rely heavily on the standardization in Bateman [1]. 


Definition 2.1 . The Gegenbauer polynomial <7*(x) is defined by 


(1 -xY-*Ci(x) 


G < A '">£d< 



(2.1) 


where 6’(A,rc) is given by 


(-l)’T(A + i)r(n + 2A) 
G ’( A ’")- 2 “n!r( 2 A)r(n + A + 1) 


( 2 . 2 ) 

□ 


Formula (2.1) is also called the Rodrigues’ formula [2, page 175]. 
Under this definition we have 


and 


0 ( 1 ) 


r(n + 2A) 
n!T(2A) 


10(01 < 0 ( 1 ). 


— 1 < r < 1 


(2.3) 


(2.4) 
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The C.egenbauer polynomials are orthogonal under their weight function (1 — x 2 ) x ~ 

/_',(! - x 2 )'-i C x k (x)C x „(x)dx = S t ,X (2.5) 

where 

h x — tt^C x ( 1) + (2 61 

n nK T(A)(n + A) l “'° j 

We will need to use heavily the asymptotics of the Gegenbauer polynomials for large n 

and A. For this we need the well-known Stirling’s formula: 


(2tt) 2 x x+i e x < r(j: + 1) < (27 t)2 x x+ ? e x > 1 

Lemma 2.2 . There exists a constant A independent of A and n such that 

A2 


(2.7) 


A ~ 1 A2 


C x (l)<h x <A- 


C x ( 1) 


(« + A) ~ n v 7 ~ n ~ * (n + A) 
The proof follows from (2.6) and the Stirling’s formula (2.7). 


( 2 . 8 ) 


□ 


We also need the following lemma, which is easily obtained from the Rodrigues’ formula 

( 2 . 1 ): 

Lemma 2.3 . For any A > 1 we have: 

<7(A, n) 


A 

dx 


- [o - 


(2.9) 


G’(A — 1, n + 1) 

The proof follows from taking one derivative on both sides of the Rodrigues’ formula 
(2.1), and then using it again on the right hand side. 


Finally, we would need to use the following formula [2, page 176]: 


□ 


<£(*) = 


1 


2(n + fi) dx 


-rJCAAA - c^{x)} 


( 2 . 10 ) 
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3 Truncation Error in a Sub- interval 


Consider an arbitrary L\ function f(x) defined in [— 1, 1]. Suppose that the first 0 < k < A 
Gegenbauer coefficients, based upon the Gegenbauer polynomials ( f. (a;) with the weight 
function (1 — for any constant (i > 0, over the full interval [—1, 1], are given: 

/-(*) = T£(l - **r *<*(*)/(»)*. 0<k<N (3.1) 

We are interested in finding the Gegenbauer expansion of f(x), with A ~ N , based on a 
sub-interval [n,fr] C [-1, 1]- We start by introducing the local variable £: 

Definition 3.1 . The local variable £ is defined by 

x = #(£) = + 5 (^*2) 

where 

f 6 = ^ (3-3) 

Z Z 

Thus when a < x < b y — 1 < £ < 1 . 

□ 


We consider functions /(x) satisfying 
Assumption 3.2. |/ A '(fc)| < A independent of k. 


□ 


We remark that if /(x) is an L\ function this assumption is fulfilled. 

Since we know the first N + 1 Gegenbauer coefficients, f^{k) for 0 < k < N, we define 
the Gegenbauer partial sum: 

/;w = E/'Wc;( x) (3.4) 

fc= 0 

Note that /£( x) docs not converge fast to f(x) if there exist discontinuities inside the domain. 


5 


The function f(x) has also a Gegenbauer expansion in a sub-interval [a, 6], with A ~ N. 
With £, e and 6 defined in (3.2)-(3.3), we have 

OO 

/(<« + *) = £ /, J (')C, A (0, - 1 < f < 1 (3.5) 

/=0 

where the Gegenbauer coefficients f c A (l) are defined by 

/ A ( o = ^ jji m 

Of course, we do not have / A (/) at our disposal, but only an approximation based on the 
Gegenbauer parial sum f^(x), thus we have 

s<(i) = j\ i - {’)*-*<?(, ()/&('( + m (3.7) 

How well do g^{l) approximate /^(/)? To answer this question we define 
Definition 3.3 , The truncation error is defined by 

77i 

T E( A, m, N, e) = max | ]T(/ A ( /) - ^(OJCftOI (3.8) 

- K - 1=0 

where / A (/) are defined by (3.6) and g^(l) are defined by (3.7). 

□ 


The truncation error is the measure of the distance between the true Gegenbauer ex- 
pansion in the interval [a, 6] and its approximation based on the Gegenbauer partial sum in 
[ 1 , 1 ]- 

We first have the following lemma: 

Lemma 3.4 . The truncation error can be estimated by 

|c?0) 


(3.9) 


OO in rnn r i 

TE(\,m,N,<)< £ l/“(?)l£ -TT^V + 

q-N+\ 1=0 n l J ~ l 

Proof: From (3.6) and (3.7) we have = . £ . • . 

/ A (')-i, A (0 = ^ Jji-W-ictMM + 6) -/?,(<( + S)W (3-io) 
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(3.11) 


Substituting (3.10) into (3.8), recalling (2.4) and 

oo 

M + i) -/»« + *)= E rWM + *) 

<y=N+l 


we obtain (3.9). 


□ 


For simplicity of notations we denote: 

f » =/'/!- ey-tctawM + m 

In order to estimate this term we start with the following 
Lemma 3.5 . If we denote 

/r-M 

rKi _ _Ll 

9 G(\,l) 

where G(\,l) is defined by (2.2), then we have the following recursive formula: 


(3.12) 


(3.13) 


_ l 

9 2(9 + li)t 


A— l,f+l tA— 1,/+1 


r A- 1,1+1 _ T 


<7 + 1 


A > 1, q > 1 


(3.14) 


Proof: By the definition of / ^ in (3. 13)-(3. 12), we have 

i y = fj l - 

= ah) ^C^[CUM + 0 - c;.M + *> 




i 


2G(A, /)(</ + fi)e J-1 

1 


2G(A — 1, / + 1)(<7 + n)e 


J\ Tt [ (1 “ * 8 ) A "*Q A (0] [cm* + S)- CM* + «)] * 

f\{ 1 - e) x -tc&(t) [CM* + *) - CM* + *)] * 


i 


2 (fl + /0< 


A— 1,/+1 rA-l,/+l 


tA-1,1+1 _ f 
I q-\ 


<7+1 


where we have used (2.10) for the second equality; integration by parts for the third equality 
(the boundary terms vanish because of the term (1 — £ 2 ) A 2 with A > l); formula (2.9) for 
the fourth equality; and the definition (3.12)-(3.13) for the last equality. 


□ 
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We can now obtain the following recursive estimate for I*' 1 : 

Lemma 3.6 , The 7 A,/ defined by (3.13) satisfies the following estimate: 




\ j 
7 


r(q + fi + 1 - j) 

< — 7— t — max 


C ,J *’ , J < min(A, q) 


(3-15) 


t 3 i [q + fi + 1) q-j<p<g+j 

Proof: We use induction on j. The estimate is clearly valid for j — 0. Assume that it is 
valid for j = j 0 < m\n(X,q) — 1, then 

P(9 + A* + 1 jo) 


/A,/ 

7 


< 


< 


< 


< 


- Ti/ n max 

eM [q + fi -f 1) q-jo<p<q+jo 


max 


r((7 + /x + i — jo) 
e»T(q + /i + 1) q-io<p<q+jo 
T(q + n+l -j 0 ) 1 


fX—jo,l+jo I 
p I 

1 


2 (p + 


IT1 clX 

c J °r(7 + n + 1) 2(q - j 0 + n)e q-ja<P<q+j 0 

r(? + /* — io) 


f /"A— io — 1 ,f+io + 1 

[VI 

y jA-jo— 14+Jo+l 


/•A— jo — l,/+j'o + l 
7 P+1 


i rA— io — 1 >/+io + l 

i P+l 


max /A-io-U+jo+i 

- 1 <P<7+jo +1 P 


c*+ir( 9 + /i + 1) 7-io-i<p<7+io+i 

where we have used (3.14) for the second inequality. All other steps are simple inequalities. 
This finishes the induction. 


□ 


From the previous lemma we can get the following estimate: 

Lemma 3,7 . For F A,/ defined in (3.12) with A < q we have the following estimate: 

rfo-A) |G(A,/)| 


tA,/ 


< A- 


Cr(,) |G(0,/ + A)| 


(3.16) 


where A grows at most as q 2fi l . 

Proof: For simplicity and without loss of generality we assume A is an integer. Since A < q, 
we can take j — A in (3.15) to arrive at 


/y 


r (q + fi + 1 — A) 

< — r— — max 


rOJ+\\ ^ P(? 

7- S — r ^7 r max 

1 t X F(q) 7-A<p<q+A 


I _ / s (UUIA I A. _ 

C A r(7 + // + l) q-\<p<q+X I P 

By the definition (3. 1 2)-(3. 1 3), we have, for q — A < p < q -f A, 


rO,/+A 

P 


|G(0,1 + A)| 


\Jji-e)-icu(K(<(+sn 


rCV+A 

v 
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< 


Q+.xlO^p U) f l /j _ t-2 
|G(0,/ + A)| 7-i 4j * 
r(/ + A) r(p + 2//) 1 r' 

- (/ + A)! p\ r(2fi) |G(0, / + A)| 7-i 

< , 4 ! 

- ' |G(0,/ + A)| 

where for the second inequality we have used (2.4) and for the third inequality we have used 
(2.3). Clearly A is a constant if p < \ and A grows at most as </ 2 " -1 if p > 5. Invoking 
(3.13) again we obtain (3.16). 

□ 


Using Stirling’s formula we can now easily get: 

Lemma 3.8 . For l < m < N and <7 > N, we have 

(rn + 2A)” l+2A 1 


pXJ 

Q 


< A- 


(2eA) A rn m q x 
where A again grows at most as (m + A)^ 2 ^ -1 . 

Proof: Starting from (3.16) and using the definition (2.2), we obtain: 


(3.17) 


f a./ 


< A 

< A 

< A 

< A 

< A 

< A 


r(, - A) 

|G(A,()| 

cr(,) 

|G(0,/+A)| 

1 

i-T 

r(A + i)r(/ + 2A) 

erfa) 

2 '/!r( 2 A)r(f + A + \ 

r(<j - A) 

r(A)r(/ + 2A)2 A 

fA r(<?) 

/!f(2A) 

r(« - a) 

r(A)r(m + 2A)2 a 


r(/ + A) 


m 


!F(2A) 


e A r(<7) 

( q _ A A e~ A (m + 2A) ffl+2A e-(’ n+2A )2 A 

( x q q e~ 

(m + 2A) m+2A 1 


A n q e -q 2A) 2A e _2A 


(2eA) A m m q x 

where we have used (2.2) in the second inequality; the monotonic.ity with respect to l in the 
fourth inequality; and the Stirling’s formula (2.7) for the fifth inequality. 


□ 
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We are now ready for the main theorem of this section: 

Theorem 3,9 . Let the truncation error be defined in (3.8). Let A = aeN and m = fltN 
with 0 < a, [3 < 1, then 

TE(aeN, /3eN, N, t) < A (3.18) 

where A grows at most as iV 1+2/x . In particular, if a = ft < then 

TE{aeN,aeN y N,t) < Aq tN (3.19) 


where 



Proof: The theorem follows from (3.9), the Assumption 3.2, (2.8), and (3.17). 


(3.20) 


□ 


4 Regularization Error and the Main Theorem 

The second part of the error, which is called the regularization error and is caused by using 
a finite Gegenbauer expansion based on a sub-interval [a, b] C [—1,1], to approximate a 
function f(x) which is assumed analytic in this sub-interval, has been studied in [6]. We will 
thus just qhote the result. 

We assume that f(x) is an analytic function on [a, b] satisfying 
Assumption 4.1. There exists constants p > 1 and C(p) such that, for every k > 0, 

<<•(,)£ (4-D 

□ 


max 

a<x<6 


d k f 


dx k 


(x) 


This is a standard assumption for analytic functions, p is the distance from [a, b] to the 
nearest singularity of f(x) in the complex plane (see for example [8]). Let us consider the 
Gegenbauer partial sum of the first rn terms for the function /(c£ + 8): 

m 

/*'■(() = £/*(')<?«) (4-2) 

(=0 
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with £, e and $ defined by (3.2) and (3.3), and the Gegenbauer coefficients based on [a,b] 
defined by 

/>(() = i f (l - + 8)<7«K (4.3) 

/ i ^ «/ — 1 

The regularization error in the maximum norm is defined by: 

m 

RE(\, m, e) = max /(e£ + 6 ) - X)/ A ( 0^(0 ^-4) 

/=0 

We have the following result for the estimation of the regularization error, when A ~ m 


[ 6 ]: 

Theorem 4.2 . Assume A = 7 m where 7 is a positive constant. If f{x) is analytic in 
[a, b] C [-1, 1] satisfying the Assumption 4.1, then the regularization error defined in (4.4) 


can be bounded by 

/?£( 7 m,m,f) < Aq m 

(4.5) 

where q is given by 

e(l + 2 7 ) l+27 

fl M 

(4.6) 

which is always less than 1 . 

p2 1+2 ' y 7' 7 (l +7) 1+ ' y 
I n particular, if 7 = 1 and ra = 

ftcN where (3 is a positive 

constant, then 

RE(f3N,fiN,t) < Aq tN 

(4.7) 

with 


(4.8) 


q V32 ,) 



□ 


We can now combine the estimates for truncation errors and regularization errors to 
obtain the following main theorem of this paper: 

Theorem 4.3 . (Removal of the Gibbs Phenomenon for the sub-interval case of Gegenbauer 
partial sum). 
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Consider a L i function f(x) on [—1, 1], which is analytic in a sub-interval [a, 6] C [—1, 1] 
and satisfies Assumption 4.1. Assume that the first /V + 1 Cegenbauer coefficients 


/'(*) = pr /ji-i’r'Cfi.m.ii, 


for /i > 0, are known. Let 0 < / < m be the Cegenbauer expansion coefficients, defined 

in (3.7), based on the sub-interval [a, 6], of the Cegenbauer partial sum f^(x) in (3.4). Then 
for A = ra = ficN with f) < we have 


max 


VI 

/W + «)-E». A (0 <?(*) 


/=0 


< -4 + <7r") 


(4.9) 


where 


<7r 


'27/3 V 


< 1, 


<7ft = 


'27c V 

w 


< 1 


and A grows at most as A^ 1+2m . 

Proof: Just combine the results of Theorems 3.9 and 4.2. 


□ 

We now give two remarks: 

Remark 4.3 . Comparing with the Legendre case in [6], we can see that the current proof is 
less sharp (missing a factor of f in the truncation error qj). The main loss in this sharpness 
is in the estimate (3.15). 


Remark 4.4 . No attempt has been made to optimize the parameters. 


□ 


□ 


5 Numerical Results 

In this section we give two numerical examples to illustrate our result. We will test Chebyshev 
series because these are used most often in practice. Notice that the Chebyshev polynomials 
are just Cegenbauer polynomials with /< = 0 module a constant: 7\.(r) = |(7°(x). 
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Example 5.1. We take the simple step function 


/(*)= 


1, if a < x < b 
0, otherwise 


(5.1) 


and assume that we know the first N + 1 Ghebyshev coefficients of f{x): 

f°(k) = — f (1 - x 2 y h Tk(x)f(x)dx, o < k < N 
TTCk J - 1 

where 

_ / 2, if k = 0 

Ck ~\ 1, if k > 1 


(5.2) 


(5.3) 


We then form the Chebyshev partial sum 

J%(x) = £ j°(k)Ux) (5.4) 

k= 0 

and then coinpute the approximate Gegenbauer expansion coefficient based on the sub- 
interval [a, 6] defined by (3.7): 

£(l) = i + 5)<;« (5-5) 

With these Gegenbauer coefficients, we can finally compute the uniformly accurate ap- 
proximation on [a, b] defined by 

771 

(5.6) 

/=0 

Numerical experiments (for various functions) seem to indicate that 


m 


= O.leiV, A = 0.2c (V 


(5.7) 


are good choices. Notice that in our proof we did not attempt to optimize these parameters. 

For consistency we will use (5.7) for both examples. 

For this special function (5.1), there is no regularization error. Hence all we see is 
the truncation error. In Fig. 1, left, we show the errors of a middle sub-interval [a, b] = 
[-0.5, 0.5], and in Fig. 1, right, we show that of a one-sided sub-interval [a, b] = [0, 1]. We 
can clearly see good convergence for both cases. 
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Fig. 2: Errors in log scale, f(x) defind by (5.1). [a, 6] = [—0.5, 0.5] (left) and [a, b] = [0, 1] 
(right). A = 0.2(N and m = 1. N = 20,40, 80, 160. 

Example 5.2. In the second example we take the the following function 


r , x f sin(cos(x)), if a<x<b 

M = { 0 , others! " 
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Again we assume that we know the first N + 1 Chebyshev coefficients of f{x) defined by 

(5-2). 

This time both truncation error and regularization error exist. We again pick two cases 
with middle as well as one-sided sub-intervals. From Fig. 3 we can see similar results as in 
the previous example, Fig. 1. 




Fig. 3: Errors in log scale, f(x) defind by (5.8). [a, b] — [—0.5, 0.5] (left) and [a, b] — [0,1] 
(right). A = 0.2 cjV and m = 1. N — 20,40,80, 160. 

These examples illustrate well the good convergence behavior of our approach. 

6 Concluding Remarks 

We have proven the exponential convergence in the maximum norm, of a reconstruction 
procedure using Gegenbauer series based on Cf'(x) with large A, for any L\ function in any 
sub-interval [a, b] in which the function is analytic, if we are given the first N expansion 
coefficients of this function over the full interval [ — 1, 1] based on Chebyshev, Legendre, or 
any other Gegenbauer polynomial basis. Numerical examples are also given. 
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